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Abstract 


We analyze site percolation on directed and undirected graphs with site-dependent open-site probabilities. We construct upper 
bounds on cluster susceptibilities, vertex connectivity functions, and the expected number of simple open cycles through a chosen 
arc; separate bounds are given on finite and infinite (di)graphs. These produce lower bounds for percolation and uniqueness 
transitions in infinite (di)graphs, and for the formation of a giant component in finite (di)graphs. The bounds are formulated in 
terms of appropriately weighted adjacency and non-backtracking (Hashimoto) matrices. It turns out to be the uniqueness criterion 
(N that is most closely associated with an asymptotically vanishing probability of forming a giant strongly-connected component on a 
■*—* large hnite (di)graph. 
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1. Introduction 


We are currently living in an age where many scientihc and 
industrial applications rapidly generate large datasets. The con- 
'O^nectivity and underlying structure of this data is of great inter- 
I est. As a result, graph theory has enjoyed a resurgence, be- 
coming a prominent tool for describing complex connections in 
a various kinds of networks: social, biological, technological! It- 
^ 0, and many others. Percolation on graphs has been used to de- 
I—I scribe internet stability||9][T0l, spread of contagious diseases lfTTl - 
[T3]| and computer viruses lfT4l : related models describe market 
crashes lfTSll and viral spread in social networks ifThUTSl . Gen¬ 
eral percolation theory methods are increasingly used in quan¬ 
tum information theorv lfT9lj23l . Percolation is also an impor- 
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a\ tant phase transition in its own right Il24lj27l and is well estab- 
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lished in physics as an approach for dealing with strong disor¬ 
der: quantum or classical transport Il28l430l . bulk properties of 
composite materials|[3Tl |32l, diluted magnetic transitions Il33l . 
or spin glass transitions ll34lj3^ . 

Recently, we suggested ll^ a lower bound on the site per¬ 
colation transition on an infinite graph. 


Pc > llp{H). 


( 1 ) 


Here, p{H) is the spectral radius of the non-backtracking (Ha¬ 
shimoto) matrix lTO H associated with the graph. This expres¬ 
sion has been proved||39l for infinite quasi-transitive graphs, a 
graph-theoretic analog of translationally-invariant systems with 
a hnite number of inequivalent sites. The bound Q is achieved 
on any inhnite quasi-transitive tree ll^ . and it also gives numer¬ 
ically exact percolation thresholds for several families of large 
random graphs, as well as for some large empirical graphs HTI . 
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In most applications of percolation theory, one encounters 
large, but hnite, graphs. The expectation is that the correspond¬ 
ing crossover retains some properties of the transition in the 
inhnite graphs, e.g., the formation of large open clusters be un¬ 
likely if the open site probability p is well below pc- However, 
the bound (Q tells nothing about the structure of the percolat¬ 
ing cluster on hnite graphs, and neither it gives an algorithm for 
computing the location of the crossover in the case of a hnite 
graph!?]. In particular, Eq. ([^l misses the mark entirely for any 
hnite tree where p{H) - 0. 

In this work, we construct several spectral and algebraic 
bounds for transitions associated with heterogeneous site perco¬ 
lation on directed and undirected graphs, both hnite and inhnite, 
and analyze the continuity of these bounds for a sequence of h- 
nite digraphs weakly convergent to an inhnite graph. Namely, 
for hnite digraphs, we construct explicit upper bounds for the 
local in-/out-/strong-cluster susceptibilities (average size of a 
cluster connected to a given site), the strong connectivity func¬ 
tion (probability that a given pair of sites belongs to the same 
strongly-connected cluster), and the expected number of sim¬ 
ple cycles passing through a given arc. We also construct some 
analogous bounds for inhnite digraphs, which result in non¬ 
trivial lower bounds for the transitions associated with diver¬ 
gent in-/out-cluster susceptibilities, emergence of inhnite in- 
/out-clusters, and the strong-cluster uniqueness transition. 

Our results imply that Eq. Q and its analogue for hetero¬ 
geneous site percolation on a general digraph give a universal 
bound for the strong-cluster uniqueness transition, below which 
a strongly connected inhnite cluster cannot be unique. Such a 
bound is continuous for an increasing sequence of subgraphs 
if the percolation problem on the limiting digraph has a hnite 
minimum return probability, the probability that any arc and its 
inverse are connected by an open non-backtracking path. Ei- 
nite minimum return probability also guarantees that below this 
bound, the strong connectivity decays exponentially with the 
distance, and the expected size of a strongly connected clus- 
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ter scales sublinearly with the number of vertices in a digraph. 
In comparison, the bound ([T]| applies only conditionally to the 
percolation transition proper, e.g., for a weakly-convergent se¬ 
quence of quasi-transitive digraphs of increasing size, where 
the number of inequivalent vertex classes remains uniformly 
bounded. 

The remainder of this paper is organized in four sections. In 
Sec.j^we dehne several matrices associated with heterogeneous 
site percolation and introduce other notations. Our main results 
are given in Sec. which contains bounds for hnite digraphs, 
and in Sec. where inhnite digraphs are discussed. Finally, in 
Sec. I^we compare effectiveness of different criteria in limit¬ 
ing the emergence of a giant component, an open cluster which 
contains a finite fraction of all vertices in the digraph. 

2. Definitions and notations 

We consider only simple directed and undirected graphs 
with no loops or multiple edges. A general digraph D - {'V, &) 
is specified by its sets of vertices (also called sites) W’ = W’(2)) 
and edges £ s £(£)). Each edge (bond) is a pair of vertices, 
(m, v) c £ which can be directed, m —» v, or undirected, m v. 
A directed edge m —> v is also called an arc from u to v; an undi¬ 
rected (symmetric) edge can be represented as a pair of mutu¬ 
ally inverted arcs, u ^ v = {u ^ v,v ^ u], A digraph with 
no undirected edges is an oriented graph. We will denote the 
set of arcs in a (di)graph D = Each vertex v e 'V 

in a digraph D is characterized by its in-degree id(v) and out- 
degree od(v), the number of arcs in to and from v, respec¬ 
tively. A digraph with no directed edges is an undirected graph 
0 = i'V, £). Eor every vertex in an undirected graph, the degree 
is the number of bonds that include v, deg(v) = id(v) = od(v). 

We say that vertex u is connected to vertex v on a digraph 
£), if there is a path from u - uqXov = U(, 

^ = {Mq —> U\, U\ —> 112, 1 —^ ^ JH(^-l). (2) 

The path is called non-backtracking if m,_i 4^ m,+i, Q < i < {, 
and self-avoiding (simple) if m, + Uj for 0 < i,] < {. The 
length of the path is the number of arcs in the set, £ -\'P\. The 
distance from m to v on D, din, v), is the minimum length of a 
path from u to v. We call path Q open if uq + ui, and closed 
otherwise. A closed path is a cycle; it can be non-backtracking 
or self-avoiding (simple). Connectivity on an undirected graph 
is a symmetric relation: we just say that vertices u and v are 
connected (or not). On a digraph, we say that vertices u and v 
are strongly connected iff u is connected to v and v is connected 
to m; u and v are weakly connected on D if they are connected 
on the undirected graph underlying D. A ray is a semi-inhnite 
simple path, characterized as in- or out-going according to the 
directionality of the constituent arcs. A strong ray is a strongly 
connected union of in- and out-going rays; it has the property 
that the intersection between the vertex sets is an inhnite set. 

A digraph D is called transitive iff for any two vertices u, v 
in'V = 'ViD) there is an automorphism (symmetry) of D map¬ 
ping u onto V. Digraph D is called quasi-transitive if there is a 
hnite set of vertices 'Vq c 'V such that any u e "V is taken into 


W'o by some automorphism of D. We say that any vertex which 
can be mapped onto a vertex mq 6 W'o is in the equivalence class 
of mq. The square lattice is an example of a transitive graph; a 
two-dimensional lattice with r inequivalent vertex classes de- 
hnes a (planar) quasi-transitive graph. 

A graph Q' - (W",£') is called a covering graph of ^ = 
('V,S) if is there is a function f : 'V' 'V, such that an 

edge iu',v') e £' is mapped to the edge (f(u'),f{v')) e £, 
with an additional property that / be invertible in the vicin¬ 
ity of each vertex, i.e., for a given vertex u' € W" and an edge 
(f(u'),v) 6 £, there must be a unique_edge {u',v') 6 £' such 
that f{v') - V. The universal cover ^ of a connected graph 
^ is a connected covering graph which has no cycles (a tree); 
it is unique, up to isomorphisms. The universal cover can be 
constructed as a graph with the vertex set formed by all distinct 
non-backtracking^paths from a fixed origin vq g with an 
edge (Vx, ^ 2 ) € £ if !P 2 = !Pi U m, m e £ is a simple extension of 
Vi. Choosing a different origin gives an isomorphic graph. The 
dehnition of a covering digraph is similar, except the mapping 
function / must preserve the directionality of the edges. The 
covering digraph 2) of a digraph D can be constructed from 
that of the underlying undirected graph by labeling the direc¬ 
tionality of the corresponding edges. 

2.7. Heterogeneous site percolation 

Consider a connected undirected graph 0. We dehne het¬ 
erogeneous site percolation on ^ where each vertex v e 'ViQ) 
has an associated probability pv,0 < Pv < i. A vertex is chosen 
to be open with probability p^, independent from other vertices. 
We are focusing on a subgraph Q0 induced by all open ver¬ 
tices on 0. Eor each vertex v, if v is open, let C(v) c g' be the 
connected component of 0' which contains the vertex v, oth¬ 
erwise C(v) - 0. If Civ) is inhnite, for some v, we say that 
percolation occurs. Denote 

fl. = 0v(^,{/^)) = P|p)(|C(v)| = cx,) (3) 

the probability that C(v) is inhnite. Clearly, for any pair of ver¬ 
tices u and V, 0^, > 0 iff > 0. 

Similarly, introduce the connectivity function, 

T„,v = T„,„(^, {p}) = P|^)(m e C(v)), (4) 

the probability that vertices u and v are in the same cluster. Eor a 
pair of vertices u, v separated by the distance diu, v), Tu,v can be 
bounded by the probability that v is in a cluster of size diu, v)-i-1. 
Thus, in the absence of percolation, t„_,, —> 0 when diu, v) —> 
00 . The reverse is not necessarily true. 

Yet another measure is the local susceptibility, 

Xv=Xvi&,{p])^^^p)i\Civ)\), ( 5 ) 

the expected cluster size connected to v. Equivalently, local 
susceptibility can be dehned as the sum of probabilities that 
individual vertices are in the same cluster as v, i.e., in terms of 
connectivities, 

Yv = ^ Tv,,„. (6) 

ue'V 
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If percolation occurs (i.e., with probability 6,, > 0, |C,,| = oo), 
then clearly = oo. The reverse is known to be true in the case 
of heterogeneous site percolation on quasi-transitive graphs l42l 
|43]| : X\' - °° can only happen inside or on the boundary of the 
percolating phase. 

An important question is the number of infinite clusters on 

in particular, whether an infinite cluster is unique. For in¬ 
finite quasi-transitive graphs, there are only three possibilities: 
(a) almost surely there are no infinite clusters; (b) there are in¬ 
finitely many infinite clusters; and (c) there is only one infi¬ 
nite cluster ll44H46l . This is not necessarily so for more general 
graphs. Notice that when the infinite cluster is unique, the con¬ 
nectivity function is bounded away from zero, t„.,, > 0„0i, > 0. 
In addition, uniqueness of the infinite cluster implies divergence 
of the local self-avoiding cycle (SAC) susceptibility TSAc(^i)j 
the expected number of distinct simple cycles passing through 
the arc a on the open subgraph £)'. In the case of homogeneous 
percolation on undirected transitive graphs, such a relation is 
given by Theorem 3.9 in Ref. Il47l . attributed to O. Schramm. 

When the open-site probabilities are equal for all sites of 
an infinite graph, - p, v e 'V (this is homogeneous site 
percolation), one defines the critical probabilities pc and pj, 
respectively associated with formation of an infinite cluster and 
divergence of site susceptibilities. There is no percolation, 9v - 
0, for p < Pc, but 0,, > 0 for p > pc. Likewise, Xv is finite for 
p < Pt but not for p > pj- A third critical probability, pu, is 
associated with the number of infinite clusters. Most generally, 
we expect px < Pc ^ Pu- For a quasi-transitive graph, one 
has BbI 

0 < Pj ^ Pc < Pu- (7) 

Here, p„ is the uniqueness threshold, such that there can be only 
one infinite cluster for p > p„, whereas for p < p„, the num¬ 
ber of infinite clusters may be zero, or infinite. For a degree-r 
regular tree Tr with r > 3, pu - I, Pc - ^lir - 1), while for 
hypercubic lattice, Z®, p„ = pc- 

2.2. Percolation on a general digraph 

There are several notions of connectivity on a digraph, and, 
similarly, there are several percolation transitions associated 
with a digraph Its). For any given configuration of open vertices 
on a digraph D (which induce the open digraph D') we intro¬ 
duce the strongly-connected cluster which includes v, Cstr(v) £ 
£)'. Similarly, one also considers an out-cluster Cout(v) £ ID' 
and an in-cluster Cin(v) £ £)', formed by all sites which can 
be reached from v moving along or opposite the arcs in 
respectively. Finally, there is also a weakly-connected cluster 
CundCv) formed on the undirected graph 0' underlying £)'. For 
each of these cluster types, we introduce the quantities analo¬ 
gous to those in Eqs. ([^, Q, and Q, e.g., the probability 0str(v) 
that V is in an infinite strongly-connected cluster, the strongly- 
connected susceptibility Tsti ('')> the two sided (strong) connec¬ 
tivity Tstr(M, v) which implies a path from m to v and one from v 
to u must both be open, and the directed connectivity from 
U to V. 


2.3. Emergence of a giant component 

In network theory, a percolation-like transition on a finite 
graph is usually associated with the emergence of a giant com¬ 
ponent, an open cluster which contains a finite fraction of all 
vertices in a graph. The transition is sharp and it is well under¬ 
stood in various ensembles of random graphs and digraphs, see, 
e.g.. Refs. Il49»54l. 

2.4. Matrices associated with a digraph 

For any heterogeneous site percolation problem on a di¬ 
graph with the adjacency matrix A, we associate the following 
three matrices: weighted adjacency matrix 

-i 12 1/2 

[Ap]ij - Pj AijPj (no summation), (8) 

weighted line digraph adjacency matrix Lp, and weighted Ha- 
shimoto matrix Hp. 

Weighted line-digraph adjacency matrix. For any digraph D, 
the line digraph ll55l X. = X{D) is a digraph whose vertices 
correspond to the arcs of the original digraph WiX) - -Xl(D), 
and it has a directed edge (a ^ b) e -XliX) between vertices 
a - i ^ j and b - /' —> 1 iff / = f (that is, arcs a and b, taken 
in this order, form a directed path of length two). We denote the 
corresponding adjacency matrix L, and introduce the weighted 
matrix Lp, where an entry corresponding to the directed edge 
(a ^ b) e ^iX) (a - i ^ j and b - j' ^ I are arcs in £)) has 
weight Pj-. 

(Lp)ab ^ PjSjj'. (no summation) (9) 

In the homogeneous case, pj - p for all sites, and the weighted 
matrix has the simple form, Lp - pL. Notice we used the arc 
set .XI to define the line digraph; the same definition can be used 
to associate a line digraph XiS) with an undirected graph 0. 

Weighted Hashimoto matrix. Hashimoto, or non-backtracking, 
matrix H has originally been defined for counting non-back- 
tracking cycles on graphs BOl . This matrix is the adjacency ma¬ 
trix of the oriented line graph (OLG) 'H(D) associated with the 
original (dilgraph lS^ . The OLG is defined similarly to the line 
digraph, except that the edges corresponding to mutually in¬ 
verted pairs of arcs in the original digraph are dropped. We de¬ 
fine the corresponding weighted matrix by analogy with Eq. 

(Hp)ab = Pjdjjfl - Sij), (no summation) (10) 

where a - i ^ j and b - f ^ I are arcs in the original 
digraph. Again, in the homogeneous case, pj - p for all sites, 
we recover the usual Hashimoto matrix, Hp - pH. 

Notice that in the case of an infinite digraph, the objects Ap, 
Lp, and Hp are not matrices but operators acting in the appro¬ 
priate infinite-dimensional vector spaces. Eor a locally-finite 
digraph, the action of these operators is uniquely defined, re¬ 
spectively, by the local rules ([^, (|^, and ( [T0| ). Eor convenience 
we will nevertheless refer the them as “matrices”, each time 
specifying whether the graph is finite or infinite. 
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2.5. Perron-Frobenius theory 

Consider a square n x n matrix B with non-negative matrix 
elements, Bij > 0, and not necessarily symmetric. The spectral 
radius p{B) = max,<„ |/i,(B)| and the associated eigenvectors of 
such a matrix are analyzed in the PetTon-Frobenius theory of 
non-negative matrices ll57l - l59l . In particular, there is always an 
eigenvalue Tmax = p{B), and the corresponding left and right 
eigenvectors and = d^ax^L, B^r = dmax^R, can be 

chosen to have non-negative components, > 0, ^r, > 0, 
although in general one could have p(B) = 0. Further, in the 
case where B is strongly connected (as determined by a digraph 
with the adjacency matrix given by non-zero elements of B), 
the spectral radius is strictly positive, as are the components of 
^R. For such a positive vector we will consider the height 
ratio 

y(^) = max^, r(^) > 1. (11) 

‘j 


3. Finite graph bounds 

3.1. Approach 

The derivation of Eq. ([T]) in Ref. Il39]l relied on the map¬ 
ping of the percolation thresholds between the original graph 
0 and its universal cover, a tree T locally equivalent to Q, 
Pc(0) ^ Pci'T'). Our approach in this work is mostly algebraic. 
We consider the bound ([T]) as the convergence radius for the 
infinite power series ll60l . 

OO 

M = M(H) = J^p^H\ (12) 

. 5=1 

where a matrix element of H\ [//*]„,„ gives the number of non¬ 
backtracking paths starting at site i along the arc a = i ^ j, 
ending at the arc b, and visiting s - 1 intermediate sites. Thus, 
the sum Yib ^ab is an upper bound on the average number of 
sites which can be reached starting along the arc a. Respec¬ 
tively, in an infinite graph, the convergence of the series © 
implies: with probability one any given point belongs to a finite 
cluster. Unfortunately, this argument does not limit giant com¬ 
ponents on a finite graph. Indeed, matrix H is non-symmetric, 
and convergence of the series can be highly non-uniform in s 
(see, e.g.. Ref. ED), with the norm of each term exponentially 
increasing as p^||//''|| ~ p*||//||^ for s < sq - 0(m), and only 
starting to decrease for p < IjpiH) at s > sq- Thus, formal 
convergence does not guarantee the low probability of finding 
a giant component in a large but finite graph. 

On the other hand, the expansion ( [T2| i gives a convenient 
tool to study percolation. By eliminating the contribution of 
backtracking paths, and reducing over-counting, one can get 
bounds tighter than what would be possible in a similar ap¬ 
proach based on the adjacency matrix. Now, the original prob¬ 
lem of percolation on an undirected graph is substituted by the 
problem of percolation on a directed graph, the OLG whose ad¬ 
jacency matrix is the non-backtracking matrix H. The same for¬ 
malism can be used to bound percolation on digraphs, a prob¬ 
lem of high importance in network theory ifTTHTSl . 


з. 2. Spectral bounds on susceptibilities 

In the following, we only construct bounds for the out-cluster 
susceptibilities. The corresponding bounds for in-cluster sus¬ 
ceptibilities, can be obtained by considering the trans¬ 

posed matrices, and Hj,, respectively. 

Theorem 1. Consider heterogeneous site percolation on a fi¬ 
nite strongly-connected digraph D. Assume that the spectral 
radius of the weighted adjacency matrix satisfies p(Ap) < 1, 
and let ^r be the corresponding right PF vector. The out-cluster 
susceptibility for an arbitrary vertex v e 'VfiD) satisfies: 

Uout(v) < Ci{0 = max max < yiO- (13) 

l-p{Ap) uc'V vcW 

Proof. Consider the alternative definition ^ of the suscepti¬ 
bility Xv as a sum of connectivities. Any given site u is in 
the out-cluster of v iff there is an open path leading from v to 

и. The corresponding probability can be estimated using the 
union bound, with sum of probabilities over self-avoiding paths 
upper-bounded by matrix elements of powers of the matrix Ap. 
Namely, the upper bound for the susceptibility ;^fout(r) reads: 


Uout(v)<py^^ ^ [{ApYfup'J^- (14) 

s=Q ue'ViD) 


For each u, we replace plp < ^ru/ minf^Ril which re¬ 
duces powers of Ap to those of p{Ap), thus 

1/2^ 

Uout(v) < [1 -p(Ap)]-^ (15) 

mini ^Ri/pf 

The uniform bound ( [T3] l is obtained by maximizing over v. □ 

In the case of an undirected graph, or a digraph with some 
undirected edges, we can try to construct a better bound by con¬ 
sidering only non-backtracking paths, with the corresponding 
probabilities counted using the weighted Hashimoto matrix FIp. 
The argument is simplest when the OLG of the original graph 
is also strongly connected. We have 

Theorem 2. Consider heterogeneous site percolation on a fi¬ 
nite digraph ID with a strongly connected OLG. Assume that 
the spectral radius of the weighted Hashimoto matrix satisfies 
p(Hp) < 1, and let rjR be the corresponding right PF vector. 
The local out cluster susceptibility satisfies: 


Uout(r) ^ Pmax + 


ClipR) 

^-piHpY 


Cifi) = 


maX(i Pn fj^i Tjn^i 

minb^u'^jPblPj 


(16) 


The proof is analogous to that of Theorem[D except that the 
sum ( [T^ is replaced by a similar sum in terms of the weighted 
Hashimoto matrix, with an additional summation over all arcs 
leaving a chosen vertex v. We see that C 2 (? 7 ) < odjnaxT(?7), 
where od^ax is the maximum out degree of a vertex on D. 

One possible set of sufficient conditions for an OLG to be 
strongly connected is given by the following Lemma: 
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Lemma 3. Consider a strongly-connected digraph T). The 
corresponding OLG 'HitD) is also strongly connected if either 
of the following is true: (a) D has no undirected edges, or 
(b) D remains strongly connected after any undirected edge 
i <r^ j € &{D) is replaced by either of the two arcs, a = i ^ j 
or d = j ^ i. 

Proof To connect the arcs a — i ^ j and b - u ^ v from 
take a directed path V ^ {j ^ ji,ji j 2 ,..., jm-\ ^ 
m) from j to u. If ji - i [in the case the inverse of a is also 
in the arc set, a - j ^ i e replace the first step by a 

directed path from j to i which does not include the arc j 
i. If needed, do the same at the other end, and remove any 
backtracking portions in the resulting directed path. □ 


In the special case of a symmetric digraph D correspond¬ 
ing to an undirected graph 0, the condition (b) in Lemma is 
equivalent to ^ having no leaves llWI (degree-one vertices). 

For completeness, we also establish the relation between 
the spectral radii of the matrices Ap and Hp (this is an extended 
version of Theorem 1 from Ref. ||39l; see also Ref. ED): 

Statement 4. (a) The spectral radii of the matrices (|^, and 
( |7QP corresponding to heterogeneous site percolation on a finite 
strongly-connected digraph D satisfy p(Ap) - p(Lp) > p(Hp). 
(b) If D has no undirected edges, all three spectral radii are 
equal, (c) If the OLG of a finite digraph D is strongly connected 
andpiHp) — p(_Ap), then D has no undirected edges. 


Proof. Consider a right PF eigenvector of the matrix Hp with 
non-negative components, qa ^ 0, where a = i ^ j goes 
over all arcs in the digraph £). Such an eigenvector can always 
be constructed, whether or not the corresponding graph, OLG 
of £), is strongly connected ll57l459l . According to the defini¬ 
tion ( [T0| ), the corresponding equations are 

1: j->leJi{D) 


where A = p{Hp) > 0. The last term in Eq. ( [T7] l removes the 
backtracking contribution; it is present iff the edge (/, j) 6 £(£)) 
is an undirected edge (i j). Denote xj = ^fp~j Vhi 
and y, = -fpl P- 0’ where in the latter case we are 
only summing over the neighbors j connected to i by undirected 
edges; y, = 0 if there are no undirected edges incident on i. In 
Eq. (17 1 , fix i, sum over j, and multiply by -fpi', this gives 


Axi = 


j 


(18) 


or just Ax - ApX - y. Denote f a left PE vector of Ap, such 
that f^Ap - p{Ap)f^. By assumption, the original digraph is 
strongly connected, and all p, > 0; therefore vector f is unique, 
up to a normalization factor, which can be chosen to make com¬ 
ponents of f all positive. Multiply the derived Eq. ( fTSl l by and 
sum over i e 'ViDf, this gives 

Afx = fApX - fy = p{Ap)fx - fy < p{Ap)f x. 

Notice that > 0 while Xi > 0 and vector x O', thus x > 0, 
which proves p{Ap) > A = p{Hp). The same argument for 


the line digraph adjacency matrix Lp produces identical equa¬ 
tions with y —> 0 ; this gives El p{Lp) - p{Ap) and completes 
the proof of (a). Similarly, the second term in the RHS of 
Eq. ([T§ is absent in the case where D is an oriented graph 
with no undirected edges (in this case any path is automatically 
non-backtracking); this proves (b). Further, if the OLG TT is 
strongly connected, we have 77 ^ > 0 , the components of the 
vector X are all positive, Xi > 0. In this case, p{Ap) > p{Hp) 
unless D has no undirected edges, which proves (c). □ 

Since Theorem was derived by counting only non-back- 
tracking paths, one could expect the corresponding bound to be 
tighter whenever some undirected edges are present in £(£)). 
Indeed, in this case p{Hp) < p{Ap), and, assuming p{Ap) < 1, 
the denominator in Eq. ( [Th] ) is larger than that in Eq. ( [T3] l. How¬ 
ever, for sufficiently small p(Ap), the relation between the two 
bounds is determined by the coefficients Ciifu) and C 2 iqR), and 
these depend strongly on the graph. In particular, in the case 
of homogeneous percolation with on-site probability p ona d- 
regular undirected graph, Ci(^r) = p while C 2 (riR) - dp^, so 
that, indeed, the second bound is tighter. On the other hand, 
02 ( 77 ) can be infinite when the OLG is not strongly connected— 
e.g., when the graph contains a degree-one vertex. [Notice that 
the structure of the coefficients Ci(^) and 02 ( 77 ) ensure that they 
remain bounded when one of the on-site probabilities p,, be¬ 
comes very small, as long as the graph remains well connected 
after the removal of the vertex v and incident edges.] 

Einally, in the case of a strongly-connected oriented graph 
with no undirected edges, where Statement|^gives A = p{Ap) - 
p{Hp), one gets 02 ( 77 ^) = ACffR), so that the bounds in Theo- 
rems[2and|^become almost identical, as expected. 

In the following subsection, we give several bounds on ver¬ 
tex connectivities. An exponential decay of connectivity im¬ 
plies a sublinear scaling of the expected size of an open cluster 
with 7 ? s 1^1, and, by the Markov’s inequality, an asymptot¬ 
ically zero probability of a giant open component which con¬ 
tains cn or more vertices, with any c > 0. This can be seen 
from the following 


Lemma 5. Consider heterogeneous percolation on a digraph 
D with n — \'V(fD)\ vertices. For a chosen site v, assume 
that the directed connectivity decays exponentially with the dis¬ 
tance, i.e., there are some constants C > 0 and p < 1 such that 
T,,_„ < Cp"'*''’“\ while the number of sites at distance mfrom v 
is exponentially bounded, |(m € 'V : d(v,u) - 777 )! < C'A’”, for 
some C > 0 and A > 1. Then, for any q such that x = Apt < 1, 


L'out(v) <77* 


CiCfIt 


(19) 


Proof. Consider Eq. as a 1-norm of an 7?-component vec¬ 
tor of connectivities. Eq. ( [T^ follows immediately from the 
Holder’s inequality, ||£7||i < ||£7||^7?^“'^^, where the ^-norm is up¬ 
per bounded in terms of an infinite geometrical series with the 
common ratio x. □ 


Note; for a digraph whose underlying graph has a maximum 
degree limax, A < i7max - 1- 
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3.3. Spectral bounds on connectivity and SAC susceptibility 
It is easy to check that whenever the height ratio of the cor¬ 
responding PF vector is bounded, all connectivity functions de¬ 
cay exponentially with the distance if p{Ap) < 1. However, 
since only two vertices are involved, we can get a more general 
statement; 


where jo = i, and j\ - je = j; we assume V c Denote 

t 

P = Pi'P)^ A'^YYpj. and xi = ^ rp^i. (25) 

i=l ti^leJi(X>) 

Then the following inequalities are true: 


Theorem 6. Consider heterogeneous site percolation on a fi¬ 
nite digraph D characterized by the weighted adjacency matrix 
Ap with the spectral radius p = p{Ap) < 1. Then, for any pair 
of vertices, the directed connectivity either from u to v or from 
V to u is exponentially bounded, 

Tu,v < (1 - or T,, „ < (1 - (20) 


where d(u, v) is the directed distance from u to v. 


Proof. The statement of the theorem follows immediately if u 
or V are not strongly connected, as in this case either or both 
directed connectivities are zero. When u and v are strongly con¬ 
nected, the corresponding components of the right PF vector of 
Ap are both non-zero, fu > 0, fv > 0. The argument similar to 
that in the proof of Theorem[2gives the bound 


Tu,v < 


^P^ 

1 - p ' 


( 21 ) 


The pair with smaller prefactor proves (201. 


□ 


Theoremj^can be interpreted as a bound on the strong con¬ 
nectivity, and an associated bound on strong-cluster suscepti¬ 
bility from Lemma 


Corollary 7. Consider heterogeneous site percolation on a fi¬ 
nite digraph ID with n = fV{D)\ vertices, characterized by the 
weighted adjacency matrix Ap with the spectral radius p = 
p(Ap) < 1. Then, for any pair of vertices, the probability that 
they are in the same strongly connected cluster is exponentially 
bounded. 


XjPjP 

“ Pj + ^P' 


(26) 

(27) 


Proof. First, notice that q satisfies the eigenvalue equation ( [T7] i. 
Using it repeatedly along V, write 


Ph 


dh~>ji - ^ Pji^ji 


PhPji 


A2 




Ph 


Af 




(28) 

which gives Eq. ((26|. Second, notice that the sum in the RHS 
of Eq. (17 1 equals xj. Combine Eq. (17 1 with a similar equation 
for the inverted arc j —» i which contains x,. Together, the two 
equations give (we assume A^ + PiPf) 

Axj - piXi Axi - pjXj 

" Pi ^2. „ „ ’ ^i^< = Pi 


• PiPj 


d2- 


• PiPj 


Substite in Eq. (26 1 to obtain 


PiXi 


A + Ppi PjXj 


(30) 


'P'- PiPj Pj + AP X^- pipj ’ 
where the denominator is preserved for its sign. The lower 


bound (27 1 is obtained by substituting back into 


□ 


Eor every vertex j of a digraph D with a strongly-connected 
OLG, let us define the minimal return probability, 

Pj = min I —, min max P(!P) \, a e (31) 


Tst,(M, V) < (1 - d(v,u)] (22) 


In addition, if the underlying graph ofD has a maximum degree 
t/max. Strong-cluster susceptibility satisfies. 


Xstx{v)<Cn'' 


C = 


1 

1 -p 



(23) 


for some q > 0 such that x = (r/max ~ l)p^ < 1- 

We now would like to prove a tighter version of Theorem 
l^which involves non-backtracking paths and the weighted Ha- 
shimoto matrix Hp, see Eq. (101. To this end, we first quantify 
the strong connectedness of the OLG of the digraph D: 


Lemma 8. Let q be the right PF eigenvector of the matrix Hp 
corresponding to the positive eigenvalue A = p(Hp) > 0. Con¬ 
sider a pair of mutually inverted arcs a = i ^ j and a = j ^ i, 
[a, d} C IR{D), connected by a length-I simple path. 


P = {jo ^ ji> ji ji,---, je-i je, je jo), (24) 


where the minimum is taken over all arcs leading to j, V is 
a non-backtracking path (24i connecting a = I ^ j and its 
inverse, a = j ^ I, and P((P) is defined in Eq. (25 1 . Eor a - 
i j such that the inverted arc does not exist, a i IR{D), we 
should use P - PjjA, see Eq. (17 1 . Notice that thus defined Pj 
is a local quantity; we expect it to be bounded away from zero 
for (di)graphs with many short cycles. 

We can now prove 


Theorem 9. Consider heterogeneous site percolation on a fi¬ 
nite digraph D characterized by the weighted Hashimoto ma¬ 
trix Hp with the spectral radius p = p(Hp) such that 0 < p < 1. 
Let the OLG of D be also strongly connected, with a non-zero 
minimal return probability, Pmin = min Pj > 0. Then, for 
any pair of vertices, the directed connectivity either from i to j 
or from j to i is exponentially bounded. 


Hi < 


rf(ij)-l d(j,i)-l 

^ (2^™!n)> Hi < ^ 


(1-p) 


(1-p) 


(2^;!„)- (32) 
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Proof. Start with the union bound for the directed connectivity, 
formulated in terms of non-backtracking paths from i to j, 

Z Z Z ^33) 

m>d{ij)—l l:a=i^IeJ{ V\b=l'^jeJi 


where we assume i 4^ j, so that d(i,j) > 1. Construct a further 
upper bound by introducing the factor [cf. Eq. @1 


JP 


> 1 , 


= min ni-. 
'nun " 


XJPjpj 

Pj+pp/ 


(34) 


and extending the summation over b to all arcs. This replaces 
the matrix element [H'”]ab with p"', 

nj<p,pj Z Z ^35) 


The summation over arcs a from i gives Xi in the numerator, see 
Eq. ( [25] l, while the lower bound ( [34) l gives Xj in the denomina¬ 
tor. We obtain 


Ti., 



The uniform connectivity bound ( [3^ is obtained after the sum¬ 
mation using the inequalities pjjp > Pj > P^m^ Pi ^ 1- □ 

This gives a conditionally stronger version of Corollary]^ 

Corollary 10. Consider heterogeneous site percolation on a 
finite digraph D characterized by the weighted Hashimoto ma¬ 
trix Hp with the spectral radius p = p(Hp), 0 < p < 1. Then, 
if D satisfies the conditions ofLemma^ the probability that a 
pair of vertices are in the same strongly-connected open cluster 
is exponentially bounded, 

„min[£/(M,v’), d{v,u)]-l 

Tsifu, V) < - (2P-Jj. ( 36 ) 

Exponential decay of connectivity with the distance also 
implies a sublinear scaling of the expected size of a strongly- 
connected cluster (Lemma 1^: 

Corollary 11. Consider heterogeneous site percolation on a di¬ 
graph D with n = YV(fD)\ vertices, characterized by the weight¬ 
ed Hashimoto matrix Hp with the spectral radius p = p(Hp), 
0 < p < 1. Let ihe maximum degree of the undi¬ 

rected graph underlying D, and q > I satisfy the inequality 
X = (t/max ~ l)p^ < 1- Then, if the graph satisfies the condi¬ 
tions of Lemma we have the following bound for strongly- 
connected cluster susceptibility. 




(2P-i ) 

' min'' 


(1 -p)(l -x)'/?' 


(37) 


This bound guarantees that in a large digraph, a giant strongly 
connected component occurs with asymptotically zero proba¬ 
bility as long as the corresponding prefactor remains bounded. 

We conclude this subsection with the following universal 
bound on the expected number of SACs through a given arc; 


Theorem 12. Let Hp be the weighted Hashimoto matrix {10\ 
for heterogeneous site percolation on a finite digraph D. Then, 
if the spectral radius p(Hp) < 1, the SAC susceptibility for any 
arc a £ J{(fD) is bounded. 


TsAc(a) < [1 -p{Hp)Y 


(38) 


Proof We hrst consider the case of a digraph whose OLG is 
strongly connected. The expected number of SACs, Ns{a), of 
length s passing through the arc a on £)' is bounded by 


Ns(a) < [Hp]aa (no summation.) (39) 

Let ?7 be a right PE vector of Hp with all-positive components; 
qa > a e The matrix elements of Hp being non-negative, 
the RHS of Eq. ( (39l l can be further bounded by 

= W7/p)]^ (40) 

where there is no implicit summation over the chosen a e ^(D). 
Summation over s gives Eq. p8] l. 

Eor a general digraph £), not necessarily connected, we no¬ 
tice that the bound p8] ) is independent of the actual values of the 
components of the PE vector q. A general hnite digraph can be 
made strongly connected by introducing additional vertices and 
additional edges connecting these vertices to different strongly 
connected components of £), with arbitrarily small probabilities 
Pi - e > 0 for these vertices to be open. The statement of the 
Theorem is obtained in the limit e ^ 0. □ 


3.4. Matrix norm bounds 

Eor digraphs where the height ratios Jr and ji may be large, 
we give a weaker general bound for the local in-/out-cluster 
susceptibilities; 

Theorem 13. Let the induced one-norm of the weighted Hashi¬ 
moto matrix corresponding to heterogeneous site percolation 
on a digraph D satisfy ||//p||i < 1. Then, the in-cluster suscep¬ 
tibility for vertex j is bounded, 

;rin(;)< l+id(;)(l-||i/;,||i)-', (41) 


where id( j) is the in-degree of j. 

We note that ||//p||i equals to the maximum column weight 
of Hp. The analogous result for the out-cluster susceptibility is 
obtained by considering the transposed matrix H^, which gives 
the maximum row weight of Hp (also, ||//^||i = ||77|U). 

Proof of 77ieorew[73] We start with a version of Eq. ( [T4| ) for 
in-cluster susceptibility. 


Xmii) < Pj + 


Pi'^YHpYYuPj 

s=0 


(42) 


where the hrst term accounts for the starting point j, summation 
over the arcs v s / —> /' and u = / —> j is assumed, and the 
expression inside the norm is a vector with non-negative com¬ 
ponents labeled by the site index i. Statement of the Theorem 
is obtained with the help of the standard norm expansion, also 
using Pi < 1. □ 
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Notice that Theorem[T^works for finite or infinite digraphs. 
This implies that in- and out-cluster site-uniform percolation 
threshold in an infinite (di)graph satisfy, respectively, 

> ii//iir‘, (43) 

pioui) > > ii//^ii7'= ||//||-'. (44) 

Thus, Theorem [T^ gives a direct generalization of the well- 
known maximum-degree bound 16^ . 

Pc>(d^..-ir\ (45) 

to heterogeneous site percolation on a digraph. 

Notice that any induced matrix norm satisfies the inequality 
\\Hp\\ > p{Hp), so that these bounds are generally weaker than 
Eq. Q. Nevertheless, Example [^below shows that the bounds 
( [4^ and ( |44| are tight: the finiteness of the height ratio jr is an 
essential condition in Theorem[T]and Theorem|2] 

Example 1. Consider a family of girth-L random strongly con¬ 
nected oriented graphs with 21? vertices, parametrized by injout 
degrees d\ > d 2 > I [see Fig. The graphs are constructed 
from 2L directed cycles of length L, by adding D, randomly 
placed arcs from each site of the i th cycle to those of the cycle 
(i H-1), with the constraint that the in-degrees on the latter cycle 
are all the same (and equal to D, -t- 1). We set D, — d\ — I for 
i — 0,... ,L — and D, — d 2 — I for i - L,... ,2L — with the 
sites of the last cycle connecting to those of the first. 

In the degree-regular case d\ - d 2 — d, the PF vectors of 
the Hashimoto matrix have height ratios Jl — Jr - 1. which 
implies that the in-/out-cluster susceptibilities for this family of 
digraphs is bounded for p < p{Hy^ — l/d. 

More generally, for this family of digraphs, p(A) - 1 equals 
the geometrical mean of the parameters D,. In the case d\ > 
d 2 , we get p{A) - p{FI) = 1 H- [(d\ - l)(d 2 - 1)]*^^, with the 
height ratios Jr and Jr divergent exponentially with L. In this 
case the more general weaker bounds \43\ and ( |44p apply, with 
ll^lli = Halloo - dy > c/ 2 - Numerically, for c/i = 3 and c /2 = 2, 
we find = 0.346 ± 0.01, see E/g. 

The bound in Theorem is local and independent of the 
size of the graph. The following gives bounds for the in-/out- 
cluster susceptibilities averaged over all sites: 

Theorem 14. Let the q-nonn (q > 1) of the weighted adjacency 
matrix (j^/c^r heterogeneous site percolation on a digraph with 
n vertices satisfy \\Ap\\q < 1. Then the q-th power average of the 
out-cluster susceptibilities satisfies 

(46) 

^ i=l 

Proof. Write a version of the bound ( [T^ as 

00 

Tout(v) < f^lA^pUeu, (47) 

s=0 

where we introduce the vector with all-one components, £« = 1, 
u - I,... ,n. Consider this expression as an element of a vector 



Figure 2: Averages of the largest out-cluster size as a function of on-site prob¬ 
ability p for digraphs in Example^ using di = 3 and d 2 = 2. Shown are the 
results for two digraphs, with L = 15 and L = 100, see Fig. □ for a smaller 
sample. Each point is an average over 120 realizations of the open vertex con¬ 
figurations with given site probability; the corresponding standard errors are 
smaller than the point size. Lines show the linear fits to the numerical data, 
taken in the regions where the cluster size fluctuations are large (not shown). 
Inset shows the finite size scaling of the corresponding intercepts, giving the 
out-cluster percolation transition at = 0.346 ± 0.01. 


Xof susceptibilities. Taking the ^-norm of this vector and using 
the standard norm expansion, we get, in self-evident notations 


WWq = ^ f^iApYe 


q 

Vq 


Eq. ([46]l immediately follows. 


(48) 

□ 


We note that one can also derive a version of Eq. (46 1 with 
the non-backtracking matrix. However, in the cases where the 
^-norm of a large matrix can be computed efficiently, q - I, 
q - 2, and q - oo, the resulting bounds are superseded by 
Theorem 13 it is easy to check that l^ \\IIp \\2 - ||7/p||i, and 
mih = iiH^iico. 

In the special case of heterogeneous site percolation on an 
undirected graph Q where the matrix Ap is symmetric, Ap - A^, 
the two-norm of Ap equals its spectral radius, ||Ap ||2 = piApf, in 
this case we obtain a stronger version of Corollary for undi¬ 
rected graphs: 


Corollary 15. Let the spectral radius of the weighed adjacency 
matrix (|^ for heterogeneous site percolation on an undirected 
graph with n vertices satisfy p(Ap) < 1. Then cluster suscepti¬ 
bility at any site i is bounded. 


X(i) < 


l-p(Ap)' 


(49) 


This bound guarantees the absence of a giant component in 
a sufficiently large graph with p{Ap) < 1. The scaling with n 
is not an artifact of the approximation; this is illustrated by the 
following 
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Figure 1: Sample of the two-region digraph from Example[^ The digraph is constructed from 2L directed cycles of length L, by adding D/ arcs from each vertex of 
/th to randomly selected vertices of the (/ + l)st cycle, so that the in-degrees on cycle i + 1 all equal (D/ 4- 1). For 0 < i < L, choose Dj = di - for L < i <2L (last 
ring connects to first), Dj = d 2 - 1. Shown is an actual sample with L = 4, di = 3, d 2 = 2. For such graphs, at large L, in-/out-cluster percolation is determined by 
that in the denser region, ^ (assuming dy > d 2 ), see text. 


Example 2 (Percolation on a rooted tree). Consider an undi¬ 
rected r-generation D-ary rooted tree with 


n = 1 + D + ... + D'' 


- 1 


D- 1 


vertices. Homogeneous site percolation with probability p gives 
the susceptibility at the root of the tree 


{pDY^^ - 1 

;^(0) ^p[\+pD + ... + {pDf] = ly ’ ^ ■ 

pD - 1 

In the limit r ^ oo percolation threshold is pj — I ID, same as 
Pc for the infinite tree T^, d — D + 1. On the other hand, the 
spectral radius of the adjacency matrix is 

p(A) - 2 VDcos(7r/(r + 1)) < 2 Vd. 

For p — llp{A) > (in the case D > 4), at large n, ;k^(0) = 
in agreement with the bound ^^ 

4. Bounds for infinite digraphs and digraph sequences 

4.1. Bounds for infinite digraphs 

In the majority of graph theory applications, one is inter¬ 
ested in perhaps large but nevertheless finite graphs, consider¬ 
ing the percolation transition as an idealization of a finite-graph 
crossover which gets sharper as n ^ cm. It would be conceptu¬ 
ally easiest to define an infinite limit by considering a sequence 
of increasing induced subgraphs, A'iQi) c A'(Qi+i), such that 
every vertex in the infinite graph Q is eventually covered. How¬ 
ever, this is not necessarily a good idea. In addition to regular 
D-dimensional lattices (e.g., iP), or graphs which can be em¬ 
bedded in TP with finite distortions, there are many important 
graph families which are, in effect, infinite-dimensional. Such 
graphs are characterized by a non-zero Cheeger constant. 


such a graph, a boundary of a subgraph may dramatically af¬ 
fect its properties, including the nature and the location of the 
percolation transition in the infinite-size limit. 

An important property of a non-amenable graph Q is that 
even if the analog f of the PF eigenvector of the graph’s ad¬ 
jacency matrix A can be constructed, Af - (e.g., in the 

case of a quasi-transitive graph), f cannot be approximated by 
a sequence of normalizable vectors. Moreover, if one consid¬ 
ers A as a bounded operator acting in the Hilbert space P'i’V) of 
finite-two-norm vectors in the vertex space the PF eigen¬ 

value dmax lies outside of the spectrum cr(A). Indeed, if we take 
any finite subset AV c 'V, the corresponding finite-support ap¬ 
proximation p of f obtained by setting for v 6 AV, and 

zero otherwise, violates the bonds along the boundary dW. As 
a result, the ratio 

IIAf - d,naxf lb 

Wf'Wl 


does not become arbitrarily small, no matter how large the sup¬ 
port \AV\ is. In comparison, the ratio (511 can be made arbitrar¬ 
ily small in an amenable graph with hg - 0. To distinguish from 
the usual PF spectral radius (with a non-normalizable eigenvec¬ 
tor), we will denote the spectral radius of A treated as an oper¬ 
ator in l^(A^) by pp(A) = sup{|/l| : A e cr(A)), with the spectrum 


(t(A) s {A : (AI - A) * is not bounded in flA^)}. (52) 


Regardless of these complications, extending our percola¬ 
tion bounds to infinite digraphs is conceptually simple. Indeed, 
e.g., in the case of uniform undirected percolation, the key ob¬ 
servation!^ is that percolation transition pc on any graph 0 
cannot be lower than that on the universal cover pc > pc. 
Universal cover is a tree, and the percolation transition on a 
locally finite tree T can be found in terms of its branching 
number 1641 . pc - l/br(7“). It turn, the branching number can 
be bounded with the growth of the tree l^ . 


\dG'\ 

mm —^-, 

g'c.g mind^'l, \Q \ Q'\) 


(50) 


where \dQ'\ is the size of the boundary, number of edges con¬ 
necting the subgraph Q' and its complement, Q \ Q', and, e.g., 
\Q'\ is the number of vertices in the subgraph. An infinite graph 
with a non-zero Cheeger constant is called non-amenable. For 


grr = indd > 0 ; liminf d-"j5„(0l = O}, (53) 

where |5,„(0I is the number of sites at distance m from the cho¬ 
sen origin /; this definition is independent on tl^ choice of the 
initial site i. Notice that when the tree 7“ = ^ is the univer¬ 
sal cover of a graph Q, |5,„(/)| is a number of length-m non¬ 
backtracking paths starting at /; this number is the same on a 
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graph and its any cover and can be counted with the help of the 
Hashimoto matrix. 

To construct a similar bound for the case of heterogeneous 
site percolation on a locally-finite digraph, we similarly asso¬ 
ciate growth with any non-negative matrix, 

gr = max inf{d > 0 : lim inf V = o}. (54) 

b 

It is a simple exercise to check that on a strongly-connected 
digraph the infinum is independent of the chosen arc a. We 
prove the following generalization of the bound Q: 

Theorem 16. Consider heterogeneous site percolation on a 
strongly connected locally-finite digraph Q. Let on-site prob¬ 
abilities be such that the growth of the corresponding weighted 
Hashimoto matrix satisfy gvHp < 1. Then, with probability 
one, any out cluster on the open subgraph is finite. 

The proof is a variant of the first half of the proof of Theo¬ 
rem 6.2 in Ref. 1641 . 

Proof. Suppose an open out-cluster for a chosen vertex v € ^ is 
infinite. Then, there should be a simple path from v to infinity, 
entirely constructed from open vertices. Then there is also a 
directed SAW from some arc a e Ji leaving v to infinity, which 
implies that for any m > 0, the number of points reachable 
from a in m steps should be at least one. Given the probability 
flout(v) > 0 to have an infinite out-cluster starting from v, we 
have the following uniform bound for any m > 0, 

pydv max 'VlHpjab > doutiv). (55) 

By assumption gr//p < 1, thus liminf^^oo 'Zih[H'f]ab = 0 for 
any arc a e That is, there should be an increasing sequence 
m,, i - 1,2,..., such that the corresponding limit is zero. Then, 
Eq. ( [55) l gives 0out(v) = 0, QED. □ 

Similarly, to get a bound on the transition associated with 
divergent susceptibility, we need to ensure the convergence of 
the series similar to that in Eq. ( [T5| ). To this end, we introduce 
the uniformly bounded erowth %5l . 

^H = inf{d > 0 ; lim sup 4^“ sup V = O}. (56) 

m—>oo u ., 


Take e - I and gxHp)l2 < 1, and define a constant 

A = (58) 

m=0 ^ b 

The supremum is finite for any finite m by the assumption of 
a finite maximum degree. Now, use our usual bound for the 
out-cluster susceptibility, in terms of all non-backtracking paths 
from V e 'V(D), counted using the powers of Hp. Summation 
over m gives the stated bound with C = A H-/l'”»+V(l-4). □ 

In comparison, the connectivity is bounded in terms of the 
fiJl) spectral radius ppfiHp): 

Theorem 18. Consider heterogeneous site percolation on a 
locally-finite digraph D characterized by the weighted Hashi¬ 
moto matrix Hp. Then, if p = pp(Hp) < I, there exists a base 
p' < I and a constant C > (1 — p') * such that the directed con¬ 
nectivity between arcs u and v decays exponentially with the 
distance d{u, v), 

T„,, < C{p't^’^\ (59) 

Proof. Start with our usual bound in terms of the weighted Ha¬ 

shimoto matrix, 

r„,v,< (60) 

m>d(u,v) 

the series converges since p < 1 by assumption. Notice that in 
fiJl), the spectral radius of Hp can be defined as the limit, 

p = lim ||//;H0“, (61) 

where > p for every m > 0. Then, for any e > 0, there 

exists mo such that < p -h e for all m > mo. Choose 

e = (1 - p)/2, define p' = p H- e = (1 H- p)/2, and also the 
constant B — maxo<,„<mo \\Hp\\/(p')"', B > 1. The statement of 
the theorem is satisfied with C = (1 - p'fy^B. □ 

The spectral radius pp(Hp) satisfies the following bounds: 

Statement 19. Consider heterogeneous site percolation on a 
locally-finite digraph D characterized by the weighted Hashi¬ 
moto matrix Hp, and the corresponding problem on the uni¬ 
versal cover T) characterized by the matrix Hp. The following 
bounds are true: 


In particular, if // is a Hashimoto matrix of a tree, then ll64l 
Pt - 1 / ^ //. More generally, we have 


Theorem 17. Consider heterogeneous site percolation on a 
bounded-degree digraph D, characterized by the weighted Ha¬ 
shimoto matrix Hp. If ffHp < 1, then there exists a constant 
C >0 such that for any v e A^(D), ;ifout(v) < 1 -i- C od(v). 


Proof. The definition (561 implies that for any A > gi Hp and 
any e > 0, there exists mo such that for all m > mo. 


4-'" sup y [H'”U<e. (57) 

aeja(O) beJi{D) 


pp{Hp) < pp(Hp) < gr{Hp). 


(62) 


Proof. Eor any 4 4: 0 and u e IB, define the vectors with 

components ^l'”\u) = 4^'"[H”']„,, and^(M) = Ii,„>o The 

parameter 4 4: 0 is outside of the spectrum of Hp iff the series 
Tjm>o Hp I'I-'” define a bounded operator on [cf. Eq. (52 1 ]. 
Equivalently, since individual arcs form a basis of for 

any u e IB, the vector f(u) should have a finite norm ||f(M)|| 2 . 


Since the matrix elements {Hp)i,v > 0, the sum in Eq. (56 1 is the 


1-norm of b"‘fu). On the other hand, for any 4 > grHp, we can 
write 

Il^(«)ll2 < X Il^^'"’(«)ll2 ^ X (63) 


m>0 


m>0 
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where the rightmost series is convergent since it is asymptoti¬ 
cally majored by the sum of Hp!Xf'-, this proves the upper 
bound. To prove the lower bound, consider a similarly defined 
vector ^(m) on the universal cover; components are the 
sums of the non-negative components of |^(m) corresponding 
to different non-backtracking paths from u to v; we thus have 
||^(m )||2 < II^(m)I| 2 - Thus, a point outside the spectrum cr{Hp) is 
also outside the spectrum cr{Hp). □ 

Further, in the case of homogeneous percolation, the spec¬ 
tral radius pp{H) - is exactly the tree’s point spectral 

radius 1641 , where the branching number hr £) < gr // < gri/. 
Also, for any quasi-transitive digraph £), we have hr £) = gr //; 
the corresponding value can be found as the spectral radius 
p(H') for a finite graph. 

The following Theorem is a generalization of Theorem 3.9 
from Ref. HtI . attributed to O. Schramm. 


Theorem 20. Consider heterogeneous site percolation on an 
infinite locally finite digraph D, with site probabilities bounded 
from above, p^, < Pmax < 1, v e 'T’(£)). Then, if strong per¬ 
colation occurs, and a strongly connected infinite open clus¬ 
ter is unique with probability one, the SAC susceptibility is un¬ 
bounded, 

supAfsAc(a) = (64) 

aej{ 


Proof. First, uniqueness of the infinite strongly-connected open 
cluster K c !D' implies that, with probability one, K is one- 
ended; it cannot be separated into two or more strongly con¬ 
nected infinite components by removing any finite set of ver¬ 
tices. Indeed, otherwise, we would have a non-unique infinite 
cluster with a finite probability, which contradicts the assump¬ 
tion. Second, with probability one, K contains two disjoint 
strong rays. Indeed, let us assume that not to be the case. Then, 
according to Menger’s theorem, for any v e the open 

subgraph would have an infinite number of single-vertex cut 
sets separating v from infinity or infinity from v. This would im¬ 
ply Oshiv) - 0, counter to the assumption. The one-endedness 
of K implies that outside of any finite ball, the two strong rays 
must remain strongly connected with each other. This means 
that with probability one, for some a e Jl, the open strongly- 
connected cluster K contains an infinite number of simple cy¬ 
cles passing through a, which implies Eq. (641. □ 


Notice that an upper bound on probabilities p^, is an essen¬ 
tial condition. This eliminates the case of a digraph with a ray 
whose vertices all have p,, = 1 . 


Given an infinite graph Qi'V, &), we say that a sequence of 
graphs t — 1 , 2 ,..., (weakly) converges to Q near the origin 
Vo e 'V, if for any R > 0 there is to such that the radius-R vicin¬ 
ity of Vo on 0 for every f > fo is isomorphic to a subgraph of 
For heterogeneous site percolation, we require p^ to match 
on corresponding sites. For digraphs, we also require the bond 
directions to match (while using undirected distance to define 
the ball). In the following, when discussing a sequence of di¬ 
graphs, objects referring to the digraph are denoted with the 
corresponding superscript, e.g., the weighted Hashimoto matrix 
and its right PF vector 6 We will also use the 

same notation for the corresponding vector mapped to under 
the isomorphism map, e l^i~^), dropping any arcs not in 
and adding zeros for arcs not in 

We first compare the bounds for the transition associated 
with emergence of infinite cluster and divergent out-cluster sus¬ 
ceptibilities. The bound in Theoremj^is formulated in terms of 
the spectral radius piHp) of the Hashimoto matrix and a pref¬ 
actor C 2 (? 7 «), while bounds in Theorems [T^ and [17] are for¬ 
mulated in terms of the growth gvHp and uniformly-bounded 
growth ^Hp. A sufficient condition for continuity between 
these bounds is given by the following; 

Statement 21. Consider heterogeneous site percolation on an 
infinite digraph D, characterized by the weighted Hashimoto 
matrix Hp, and a sequence of finite digraphs with strongly 
connected OLGs, converging to D around some origin. Then, 
if the right PF vectors of H^p have uniformly bounded height 
ratios, yiq^p) < M, the following limit exists, and 

\\mp(Hf) ^ gfiHp) ^'^(Hp). (65) 

f—»oo ^ 


Proof. For each the corresponding universal cover 2)*'* is 
a quasi-transitive tree; a lift ff^ of the PF vector 77 ® is the eigen¬ 
vector of H^p with all positive components. For a given u e .71 
and t large enough, the conditions guarantee that the radius-m 
vicinity of m on £) is entirely within the subgraph of iso¬ 
morphic with that of 2 ); the same is true for the universal cov¬ 
ers. We can therefore construct the upper and lower bounds on 
the sum in Eq. (56 1 in terms of the right PE vectors q^‘\ 

-^[p{H^^)r < yiH';U < ( 66 ) 

yiriR) V 

or, using the assumed uniform bound on the height ratios, 

< y[H;u < Mip(Hf)r. ( 6 ?) 


4.2. Bounds for sequences of finite digraphs 


In Sections 3.2 and 3.3 we constructed several upper bounds 
for susceptibilities and connectivity in heterogeneous percola¬ 
tion on finite digraphs, formulated in terms of the spectral ra¬ 
dius of the corresponding weighted Hashimoto matrix Hp and 
the associated PE vectors qi, pr. In contrast, bounds in Sec¬ 
tion |4[T] are formulated directly on infinite digraphs. We would 
like to see the correspondence between these bounds for weakly 
convergent digraph sequences. 


Now, let us choose a subsequence with spectral radii converging 
to liminf,^ooP(77®)- Using only these digraphs in the upper 
bound (67 1 , the definition ([56]l implies 


gi Hp < liminfp(H®). 


The same calculation can be repeated for the lower bound, with 
a subsequence of graphs whose spectral radii converge to the 
corresponding superior limit; we get 

lim supp(H®) < gT(Hp) < MHp) < lim inf p{Hf). ( 68 ) 

r-^oo 


11 








This implies that the limit exists and satisfies Eq. ([65|. □ 


Similarly, for a finite digraph, Corollary ( [T0| ) gives strong 
connectivity exponentially decaying with the distance for p(Hp) < 
1, as long as the corresponding OLG is locally strongly con¬ 
nected. If we consider a sequence of such finite digraphs con¬ 
verging to an infinite digraph around some origin, we expect ex¬ 
ponential decay of strong connectivity for lim inf,^oo < 


1, with a bounded prefactor. On the other hand. Theorem 18 


gives exponential decay of directed connectivity with the dis¬ 
tance on an infinite digraph with pp{Hp) < 1, without an ex¬ 
plicit bound on the prefactor. The following gives partial corre¬ 
spondence between these results for undirected graphs; 


Statement 22. Consider heterogeneous site percolation on a 
locally-finite infinite connected graph Q, characterized by the 
weighted Hashimoto matrix Hp, and a sequence of finite graphs 
0^*^ converging to Q around some origin. If the ratios of the 
components of the right PF vectors of H^p corresponding 
to each arc and its inverse are uniformly bounded by a fixed 
M>\, 


any increasing sequence of finite subgraphs 0^'^ C C 0, 

convergent to 0 around some origin. Then, 

\imp{Hf)^pti{Hp). (72) 

f—»oo ^ 

Our final result is a bound on the SAC susceptibility for 
a weakly convergent sequence of finite digraphs, and an asso¬ 
ciated bound for the transition associated with the number of 
strongly-connected infinite clusters. We notice that the SAC 
susceptibility counts only finite-length cycles. This implies, 
most generally: 

Theorem 24. Consider heterogeneous site percolation on an 
infinite locally finite digraph D, characterized by the weighted 
Hashimoto matrix Hp, and an increasing sequence of finite sub¬ 
graphs C £)7+i) (3 ^ converging to D around some origin. 
Consider pq = lim,^oo p(i/®). Then, ifpo < 1, the SAC suscep¬ 
tibility at any arc a & Jiis bounded, 

;rsAc(«)<(l-Po)-‘- (73) 


M * < ^ < M, then liminfp(i/®) > pifiHp). (69) 

t—iOO ‘ 

’la 

Proof. Define the parity operator as in Ref. Il63l . Pat — da.t, to 
connect each arc a with its inverse a; we have PHpP - Hp. The 
condition of the theorem allows to use the components of Pq^‘^ 
as a lower or an upper bound on those of It also guarantees 
that the OLGs of the graphs 0^‘^ are strongly connected. For 
any m > 0 and v e there exists to such that for all t > to, the 
radius-m vicinity of v will be isomorphic to a subgraph of 0^‘\ 

If we introduce the vector e^ with the only non-zero component 
(equal to one) at the arc v, we can write for t > to, 

\\H;’eA\l = elPH';PH';e, = elP{HfrP{Hfre,, ( 70 ) 

where we used the same notation for the corresponding vector 
under the isomorphism map. We can now use the PF vector 
to construct the upper bound, 

WH’^e^Wl < elP{Hfr'P{Hfr'il^'^\ < M^\p{Hf)f'-, ( 71 ) 

W 

where we used the identity [Pq^'\ < M?;® twice. This shows 
that for any v e ,71 and any A > lim inf,^oo p(Hf'), the two-norm 
of Hp'evlA'” converges to zero, i.e., the operator (AI - Hpfi^ is 
bounded in thus pp(Hp) < liminf,^oop(77p'). □ 

We note that for an increasing sequence of subgraphs of 
0, p{H^X) < pp(Hp) in a non-decreasing bounded function of 
t. Also, the condition on the PF vector can be guaranteed by 
Lemma|^ This implies 

Corollary 23. Consider heterogeneous site percolation on a 
locally-finite infinite graph 0 whose OLG is locally strongly 
connected, i.e., for every arc a e ITI(0), there is a non-backtracking 
path of length at most I from a to a. Assume that the site prob¬ 
abilities are bounded from below, p,, > p^in > 0 . Consider 


Notice that this bound does not include limitations as in 
Corollary]^ Generally, for an increasing sequence, 

po = \\mp(Hf) < ppiHp). (74) 

t —»on ‘ 


It may well happen that po < pp(Hp), as in the case where D is 
a tree, cf. Example]^ 


Proof of Theorem [2?| The sequence piHf) is non-decreasing. 
By assumption, it is also bounded, thus the limit exists. For any 
s, the cycles of length s are contained in a finite vicinity of the 
original arc a; the corresponding contribution does not exceed 
Pg (see the proof of Theorem 12 1 . Summation over s gives the 
bound ( [73] ). □ 


Combined with Theorem]^ this gives the following Corol¬ 
lary (its weaker version previously appeared as a conjecture in 
Ref. IMl)- 


Corollary 25. Consider heterogeneous site percolation on an 
infinite, locally finite digraph D, characterized by the weighted 
Hashimoto matrix Hp, with on-site probabilities py < Pmax < 1- 
For an increasing sequence of finite subgraphs c C 

D weakly converging to D, let po = lim,^oo p(Hp\ Then, if 
strong-cluster percolation happens while po < 1, the percolat¬ 
ing cluster cannot be unique. 


Ifpo < 1 and yet an infinite strongly-connected cluster ex¬ 
ists on the induced subgraph £)', the system is expected to be 
below any transition associated with the number of percolating 
clusters. While such a transition is usually called “uniqueness” 
transition, we note that one-endedness of the infinite cluster is 
necessary but not sufficient to have a unique percolating cluster. 
An example could be any graph with a finite number of ends. 

Example 3. A degree-d infinite tree Td can be obtained as a 
limit around its root of the following graph sequences: (a) a se¬ 
quence of its subgraphs, t-generation trees (b) sequence of 
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graphs obtained from 7”® by pairing degree-one vertices arbi¬ 
trarily to form degree-two vertices; (c) a sequence of d-regular 
graphs obtained from by joining the leaves in pairs and re¬ 

placing any resulting pair of edges connected by a degree-two 
vertex with a single edge. In the case (a), — Qfor any t; 

in the case (b), lim ,= (d- 1)'^^ = with the 

components of the PF vectors falling exponentially away from 
the center. In the case (c), — d — I — gsiTj, and the PF 

vectors have all equal components. The sequence of subgraphs 
(a) correctly reproduces (the absence of) the uniqueness transi¬ 
tion, while the sequences of degree-regular graphs (c) give the 
percolation transition at pc - l/(t/ — 1 ). 

Example 4. Consider a (d, d)-regular locally-planar hyper¬ 
bolic graph 0d, with d identical d-sided plaquettes meeting at 
every vertex. This graph can be obtained as a limit of (a) an in¬ 
creasing sequence of radius-t subgraphs, with lim,^oo = 

pp-iH), or (b) a sequence of d-regular graphs whose spec¬ 

tral radii coincide with gr H — ^ H — d—l. The spectral radius 
of the adjacency matrix satisfies the following bounds based on 
the Cheeger constant \68\ \69V . 

2 'Jd-l < pp(A) < 2 yfd-, (75) 



Figure 3: As in Fig. I^but for the largest strongly-connected cluster size as a 
function of on-site probability p, giving the strongly-connected cluster perco¬ 
lation transition at = 0.530 ± 0.01. Each point is an average over 120 
realizations of the open vertex configurations with given site probability; the 
corresponding standard errors are smaller than the point size. For this graph, all 
other strongly-connected clusters are all of size 1 (single vertex). 


5. Discussion 


these result in 


V<i- 1 < pp(II) <1+ yfd. (76) 

For site percolation on Qd, we recover the maximum degree 
bound for percolation transition, pc — pr ^ ^l(d— 1), and get 
the bound pu> 11 'Jd — 1 for the uniqueness transition. 

Example 5. For a D-dimensional hypercubic lattice, iP, we 
have ppiH) - giH = gf// = 2D - 1 - d - 1. There is only 
one percolation transition, Pc — Pu — Pt- Our bounds for these 
transitions coincide and recover the maximum degree bound. 

Example 6. Consider strongly-connected cluster percolation 
on the oriented graphs in Example^ For any finite L, a strongly- 
connected cluster can be formed. Numerically, we get 
0.53, see Fig. p] This is a reasonable value since directed- 
cluster percolation in both regions is necessary to form a large 
strong cluster. However, the weak limit of this graph sequence is 
an infinite oriented tree. For such a tree, any strongly-connected 
cluster is limited to one site: strong-cluster percolation never 
happens. Respectively, the limiting spectral radius in Corollary 
|25|/.!f po = 0. 

Example 7. For an integer do > 2, consider site percolation 
on a graph Q constructed from a well-connected core, a large 
random do-regular graph Qo with no vertices, by connecting 
r additional otherwise disjoint edges to each vertex in 'V(^o)- 
With no large and r bounded, formation of a giant component on 
0 is governed by the corresponding transition on 0o. Spectral 
radii of the Hashimoto and the adjacency matrices of 0 are 

p(H) — do — \ and p(A) — |^c/o + (4r + c/q) ^ ]/2- We see that 
the bound in Corollary \1 5| becomes increasingly loose as r is 
increased in the region r > d^/A. 


We constructed several bounds for heterogeneous percola¬ 
tion on general graphs, directed or undirected. We obtained 
explicit expressions for finite and infinite (di)graphs, and, in 
two cases, analyzed the continuity of the bounds in the infinite- 
graph limit. Most bounds are obtained from the non-back- 
tracking path expansion, and formulated in terms of appropri¬ 
ately weighted Hashimoto matrices Hp. While in several cases 
stronger bounds (e.g., in terms of self-avoiding paths, or us¬ 
ing modified Hashimoto matrix as in Ref. ESI) may be readily 
available, one main advantage of the results presented here is 
that spectral radii and norms can be calculated efficiently. 

For a general infinite undirected graph, there are three tran¬ 
sitions usually associated with percolation; divergence of the 
cluster susceptibility, formation of an infinite cluster, and the 
uniqueness transition. Bounds for all three transitions are for¬ 
mulated in terms of the weighted Hashimoto matrix, see The¬ 
orems 17 16 and Corollary]^ These bounds also apply for 
directed- or strong-cluster percolation on digraphs. In addition, 
the condition ppfHp) < 1 guarantees that connectivity on an 
infinite digraph decays exponentially with the distance. 

In practical network theory applications, more important is 
the transition associated with the formation of a giant compo¬ 
nent. Several simple criteria for the emergence of a giant com¬ 
ponent that are commonly used in network theory rely on the 
degree distribution of a graph. First is the lower bound for 
percolation on an arbitrary graph in terms of the maximum 
degree ll62ll . see Eq. ( |45] l. While it is universally applicable, 
the issue with this inequality is that it tends to give very low 
bounds on graphs with wide degree distribution. Our Theorem 


13 gives a generalization of this bound to heterogeneous site 


percolation on arbitrary digraphs. Second is the Molloy-Reed 
criterion||3 [TOl EH which gives the percolation threshold 
on random graphs in terms of the two first moments of degree 
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distribution. It has been recently generalized to giant in-/out- 
clusters on random digraphs 1731. While these formulas are 
asymptotically exact in random graph ensembles ll52l . there is 
no guarantee: on actual networks the Molloy-Reed criterion can 
substantially overestimate or underestimate the thresholdlTl. 

The spectral radius has also been used to study percolation. 
On large dense graphs, under mild conditions, the critical prob¬ 
ability where a giant cluster emerges is very close to the inverse 
spectral radius of the adjacency matrix 1531 . Our Corollary]^ 
gives a related strict bound for emergence of a giant strongly- 
connected component in heterogeneous site percolation on an 
arbitrary digraph (a slightly stronger bound specific for undi¬ 
rected graphs is given in Corollary [T5] l. One substantial ad¬ 
vantage of these bounds is their universal applicability. On the 
other hand, bounds in terms of p{Ap) can become loose in cer¬ 
tain carefully designed networks, see Example |7] 

In comparison, a bound in terms of the spectral radius of 
the (weighted) Hashimoto matrix does not change upon the 
addition of leaves or finite trees. It is also asymptotically ex¬ 
act for tree-like graphs with few short cycles, large random 
graphs being the most important example. While such a bound 
does indeed limit the percolation transition on highly-uniform 
(e.g., quasi-transitive) (di)graphs, most generally the condition 
p{Hp) < 1 is a bound on the strong-cluster uniqueness tran¬ 
sition. On an infinite digraph, such a bound is constructed as 
the limit of the spectral radii for an increasing sequence of sub¬ 
graphs, which recovers the absence of the uniqueness transition 
on an arbitrary infinite tree. Further, for ppiHp) < 1, on an in¬ 
finite digraph, connectivity decays exponentially with the dis¬ 
tance. On a finite digraph, with p(Hp) < 1, strong connectivity 
also decays exponentially with an additional local strong con¬ 
nectivity condition needed to limit the prefactor in the bound. 
Such an exponential decay also implies sublinear scaling of the 
expected size of the largest cluster, which in turn guarantees 
that a giant strongly-connected component containing a non¬ 
zero fraction of all vertices emerges with an asymptotically zero 
probability. 
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